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Abstract— The application of the geometric method of local geoid model determination which requires the fitting 
of geometric surfaces to known geoid heights to enable geoid heights of new points to be interpolated involves 
the use of least squares technique for computation of the models' parameters. The selection of polynomial 
geometric surfaces depends on the size of the study area, the variation of the geoid heights and the number of 
measurement points. The accuracy of the geometric geoid model increases as the number of observation points 
approximates the number of geometric surface terms. But in most cases, the number of observation points is not 
considered. To this effect, this paper presents the relationship between geometric surfaces terms and observation 
points numbers and effect in the accuracy of geometric geoid models. A total of 23 points of known local 
gravimetric geoid heights were used. Two polynomial geometric (third and fifth degrees) surfaces were fitted to 
the geoid heights at various observation point numbers and compared to determine the relationship between the 
number of model terms and that of observation points and effect in the accuracy of the models. Least squares 
adjustment technique was applied to obtain the model parameters. The differences between the models and the 
known geoid heights of the points were computed and used to obtain the RMSEs as well as the accuracy of the 
models. The obtained results showed that the accuracy of the polynomial geometric geoid models tends to the 
highest as the number of measurement points approximates the number of the model terms and in a unique 
solution where the number of observation points is equal to the number of the polynomial geometric model 
terms, the model accuracy is highest. The paper recommends that the geometric method of local geoid model 
determination should be strictly applied in small areas. Where the method will be applied in considerable large 
areas, higher degrees polynomial geometric surfaces with a larger number of terms approximating the number 
of observation points should be applied. This will enable a proper fit of the polynomial surface to the known 
geoid heights, as well as high accuracy to be obtained. 

Keywords— accuracy, geometric geoid , model terms, points number , polynomial surfaces. 


I. INTRODUCTION 

The transformation of GNSS ellipsoidal heights to 
practical, orthometric heights in local areas has 
necessitated the determination of local geoid models of 
various areas. Local geoid model is determined using 
several methods such as gravimetric, gravimetric - 
geometric, astro-geodetic and geometric methods. The 
gravimetric method involves the use of either free air or 
Bouguer gravity anomalies computed from absolute 
gravity values of selected points obtained with a 
gravimeter within the study area and theoretical gravity 
obtained on the local ellipsoid adopted for geodetic 
computation in the study area as given by Eteje et al. 


(2019). The gravimetric-geometric method uses the geoid 
heights of points obtained from gravity observation and a 
geometric surface fitted to the gravimetric geoid heights 
to enable geoid heights of new points to be interpolated 
within the study area. The astro-geodetic method has to 
do with the use of astronomically and geodetically 
obtained data. The geometric method is applied in a small 
area and requires the use of levelling and GNSS acquired 
data such as orthometric and ellipsoidal heights of points 
accurately obtained with respect to reliable benchmarks 
and geodetic controls. Using the orthometric and 
ellipsoidal heights of the points, the geoid heights of the 
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points are computed (see figure 1) with (Oluyori et al. 
2018) 

N = h-H (1) 

Where, 

N = Geoid height 
h = Ellipsoidal height 
H = Orthometric height 



Fig. 1: Relationship between Orthometric , Geoid and 
Ellipsoidal Heights 
Source: Eteje et al. (2018) 


The geometric method has been applied by various 
researchers in different parts of the world. Such 
researchers include Erol and Celik (2004) which applied 
the fifth-degree geometric surface in an area of 

50 x 45km 2 in Turkey and Oluyori et al. (2018) that 
also applied the second-degree and the third-degree 
surfaces in the Federal Capital Territory, FCT, Abuja, 
Nigeria among others. 

Applying the geometric method for the 
determination of a local geoid model of an area, the geoid 
heights of selected points are normally computed from the 
orthometric and ellipsoidal heights of the points and a 
geometric geoid surface is fitted to the obtained geoid 
heights of the points to enable geoid heights of new points 
to be interpolated using the model. The selection of a 
method simply depends on the size of the study area, the 
variation of the geoid heights and the total number of 
points. Usually, two or more geometric surfaces are fitted 
to the geoid heights and their accuracy are computed and 
compared. The model with the highest accuracy, best fit 
the points as well as the study area and it is recommended 
for application in the study area. 

The fitting of the geometric surface to the geoid 
heights of chosen points in the study area requires the 
computation of the model parameters which in turn 
requires the application of least squares technique. In 
most cases, the lager the study area, the larger the number 
of points to be used for the determination of the geoid 


model. This enables proper depicting of the shape of the 
geoid model. 

It is assumed that the computed geoid heights have 
been adjusted. The application of the least-squares 
technique here is to obtain the model parameters. It is to 
be noted here that the reliability of the model depends on 
its ability to reproduce accurately the known geoid 
heights of the points. 

The accuracy with which geoid heights are obtained 
using the determined geoid model is computed by finding 
the differences between the known geoid heights of the 
points and their respective geoid heights from the model 
(model geoid heights). The differences are used to 
compute the Root Mean Square Error, RMSE as well as 
the reliability of the model. The accuracy of the model is 
usually highest when applying the least squares technique 
for the computation of the model parameters as well as 
fitting the geometric surface to the geoid heights of the 
points when the number of the chosen points is equal to 
the number of the model terms. In other words, in a 
unique solution where the number of observations is 
equal to the number of unknown parameters. Using the 
least squares technique, the differences between the 
estimates (most probable values) and observations is the 
residual and this is equal to the differences between the 
model and the known geoid heights of the points. 

Considering the fact that the highest accuracy of the 
geometric geoid model is obtained when the number of 
points is equal to the number of geometric surface terms, 
then when applying the method in a very large area which 
in turn requires large number of points, a geometric 
surface with a large number of terms should be applied. It 
is also to be noted here that, the higher the degree of the 
model, the larger the number of the model terms. But 
often time, the geometric surfaces are chosen considering 
the size of the study area and the variation of the geoid 
heights only. To this effect, this paper presents the 
relationship between geometric surfaces terms and 
observation points numbers and effect in the accuracy of 
geometric geoid models. 

1.2 Geometric Geoid Surfaces 

Geometric geoid surfaces are mathematical 
interpolation surfaces fitted to geoid heights to enable 
geoid heights of new points to be interpolated using 
variables such as geographic or rectangular coordinates of 
the points. These surfaces include plane surface, bi-linear 
surface, second-degree surface, third-degree polynomial 
and fifth-degree polynomial (Eteje et al ., 2018). The 
surface to be adopted as well as the degree and order of 
the polynomial depending on the size of the study area, 
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the variation of the geoid heights and the number of 
observation points. 

1.2.1 Polynomial Surface 

The polynomial surface used when determining 
geoid model is given by Erol and Celik (2004) and Kirici 
and Sisman (2017) as 

m k 

N (*,y) = Z H a i/y i (2) 

i=0 j=k-i 
i =0 

Where, 

a { j — polynomial coefficients 

m = degree of polynomial 
x, y = plan coordinates of point 
In applying the polynomial, the degree should be chosen 
and the polynomial should be formed for the chosen 
degree. Kirici and Sisman (2017) gave the third-degree 
polynomial surface with 10 terms as 

N — cLq + ci^X + a 2 Y + a 3 XY + a^X + 

a 5 Y 2 +a 6 X 2 Y+ a 7 XY 2 +a s X 3 + (3) 

a 9 Y 3 

The fifth-degree polynomial geometric geoid surface with 
21 terms is also given as 

N — cLq + a^X + a 2 Y + a 3 XY + a^X + 
a 5 Y 2 + a 6 X 2 Y 2 +a 7 X 2 Y + a 8 XY 2 + 
a 9 X 3 +a l0 Y 3 +a n X 3 Y+ a l2 X 3 Y 2 + (4) 

a l3 X 2 Y 3 + a u XY 3 + a l5 X 4 + a 16 Y 4 + 

a l7 X 4 Y + a n XY 4 +a l9 X 5 +a 20 Y 5 
Where, 

Y = ABS(y-y 0 ) 

X = ABS (x — x 0 ) 

y = Northing coordinate of observed station 
X = Easting coordinate of observed station 
y o = Northing coordinate of the origin (average 
of the northing coordinates) 

X Q = Easting coordinate of the origin (average of 
the easting coordinates) 

1.3 Observation Equation Method of Least Squares 
Adjustment 

The fitting of geometric geoid surface to a set of 
geoid heights requires the model parameters to be 
computed. The computation of these parameters is done 
by observation equation method of least squares 
adjustment technique. The functional relationship 


between adjusted observations and the adjusted 
parameters as given by Ono et al. (2014) is: 

L a =F(X a ) (5) 

Where, L a = adjusted observations and X a = adjusted 

parameters. Equation (5) is a linear function and the 
general observation equation model was obtained. The 
system of observation equations is presented by matrix 
notation as (Mishima and Endo, 2002 and Ono et al ., 
2018): 

V = AX-L (6) 

where, 

A = Design/Coefficient Matrix, 

X = Vector of Unknowns 
L = Observation Matrix. 

V = Residual 

The residual, V which is the difference between the 
estimate and the observation is usually useful when 
applying least squares adjustment technique for 
determination of local geometric geoid model parameters. 
Since it is equal to the difference between the model 
geoid heights and the known geoid heights of the points. 
So, it can be used as a check. 

The unknown parameter is computed as 

X=(A T A)- l A T L (7) 

where, 

(A r A) 1 = Inverse of the normal matrix 

The step by step procedures for the computation of 
geometric geoid model parameters are detailed in Eteje 
and Oduyebo (2018). 

1.3.1 Accuracy of Geometric Geoid Model 

The accuracy of a local geometric geoid model is 
obtained using the Root Mean Square Error, RMSE 
index. To evaluate the local geometric geoid model 
accuracy, the geoid heights of the points from the model 
are compared with their corresponding known geoid 
heights (geoid heights of the points to which the 
geometric surface was fitted) to obtain the residuals. The 
residual and the total number of selected points are used 
for the computation of the RMSE, as well as the accuracy 
of the geometric geoid model. The Root Mean Square 
Error, RMSE index for the computation of the geometric 
geoid model accuracy as given by Kao et al (2017), and 
Eteje and Oduyebo (2018) is 

lv T v 

RMSE = ±J - (8) 

V n 

Where, 

V = N Mcdd -N a ,„ (Residual) 
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N M dl = Model Geoid Height of Point 


N 


Known 


= Known Geoid Height of Point 


n = Number of Points 


II. METHODOLOGY 

The adopted methodology was divided into different 
stages such as data acquisition, data processing, and 
results presentation and analysis. Figure 2 shows the 
adopted methodology flow chart. 
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Fig. 2: Adopted Methodology Flow Chart 


The geoid heights used in this study were obtained by 
gravimetric means. They were part of the local 
gravimetric geoid heights obtained for the determination 
of the local geoid model of Benin City. The geoid heights 
were computed using the integration of modified Stokes' 


integral given by Featherstone and Olliver (1997). They 
were also corrected for combined topographic effect. A 
total of 23 local gravimetric geoid heights were used. 
Table 1 shows the 23 local gravimetric geoid heights used 
in this study. 


Table 1: Local Gravimetric Geoid Heights of the Points 


STASTION 

Northing 

Easting 

Free Air Geoid Height, N Corrected 
for Combined Topographic Effect 

XSU92 

257998.9800 

357763.3720 

2.086 

RR01 

257885.3227 

355124.0166 

2.420 

SR02 

253034.8393 

356093.6672 

1.978 

SR05 

245976.7564 

356615.1406 

2.802 

SR06 

244918.0916 

356628.3396 

3.266 

UU02 

265145.3515 

353468.5482 

3.498 
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UU03 

262403.8368 

354173.5295 

1.981 

UU05 

259407.1043 

355613.0973 

1.346 

UU08 

256422.9868 

355521.4167 

1.263 

AD01 

260514.8753 

359958.1194 

2.986 

AD03 

261867.2294 

361745.9231 

4.420 

AK02 

259528.7811 

356853.3277 

1.473 

AK05 

259332.1257 

362604.6963 

3.954 

MR02 

260751.5081 

356528.1658 

1.488 

MR04 

262930.8267 

360077.3193 

4.037 

MR05 

262428.2213 

361076.8116 

4.313 

SK02 

255516.1557 

357459.1723 

2.035 

SK03 

254396.4836 

358439.3812 

2.379 

SLK01 

259894.0672 

352909.3470 

0.781 

SLK03 

261813.3387 

350594.2641 

1.736 

SLK05 

264774.9356 

348869.1903 

3.357 

EK02 

257209.3523 

350068.7731 

0.983 

EK05 

252877.2407 

345740.0760 

2.516 


2.1 Data Processing 

The two polynomial geometric surfaces were fitted 
to the geoid heights at various point numbers. 
Considering the third-degree with 10 terms polynomial 
surface given in equation (3), 23 points, 22 points except 
for point SLK01, 21 points except for points SLK01 and 
UU03, and 11 points (XSU92, RR01, SR06, UU02, 


Degree 3 with 10 terms (23 Points) 




r 3.11075532680599771371 " 

a x 
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- 0.00000000000104326161 

a 7 
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a s 


- 0.00000000003086059461 

a 9 ) 


v - 0.00000000000746460180, 


UU08, AD03, AK05, MR05, SK03, SLK05 and EK05) 
given in table 1 were used. The model parameters of the 
various point numbers using the third-degree polynomial 
surface were computed with least squares technique as 
well as equation (6). The computed model parameters for 
23, 22, 21 and 11 points using the third-degree 
polynomial surface are respectively 

Degree 3 with 10 terms (22 Points) 
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Degree 3 with 10 terms (21 Points) 
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Applying the fifth-degree with 23 terms polynomial 
surface given in equation (4), 23 points, 22 points except 
for point SLK01 and 21 points except for points SLK01 
and UU03 given in table 1 were also used. The model 
parameters of the various point numbers using the fifth- 
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degree polynomial surface were as well computed with 
least squares technique as well as equation (6). The 
computed model parameters for 23, 22 and 21 points 
using the fifth degree polynomial surface are respectively: 


Degree 5 with 10 terms (22 Points) 


M 


^ 7.30603629707767289898 " 

a x 


- 0.01241421010635155914 

a 2 


-0.00921069831052644470 

a 3 


0.00000770915077418606 

a 4 


0.00001178307394182408 

a 5 


0.00000686623538125153 

a 6 


- 0.00000000000001862531 

a 1 


0.00000000247293704077 

a s 


- 0.00000000782726899965 

a 9 


- 0.00000000496385593976 

a io 

= 

-0.00000000150183843005 

a xx 


- 0.00000000000043701290 

a 12 


0.00000000000000007432 

a n 


- 0.00000000000000007645 

a X4 


0.00000000000176869774 

a \5 


0.00000000000077479601 

a i6 


0.00000000000009718472 

a xl 


0.00000000000000000336 

a n 


- 0.00000000000000010210 

a 19 


- 0.00000000000000003909 

V^20 ) 


v - 0.00000000000000000077, 


Page | 1186 


www.ijeab.com 





















International Journal of Environment, Agriculture and Biotechnology (IJEAB) 
https://dx. doi. org/10.22161/ijeab.4444 


Vol-4, Issue-4, Jul-Aug- 2019 
ISSN: 2456-1878 


Degree 5 with 10 terms (22 Points) 
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Fig. 3: Contour Map of Known 
Geoid Heights Using 23 Points 


The reliability, as well as the accuracy of the two 
polynomial geometric surfaces using the various point 
numbers, were computed by finding the differences 
between the models' geoid heights of the points and their 
corresponding known (gravimetric) geoid heights. The 
computed differences and the total number of points for 
each case were used to compute the Root Mean Square 
Error, RMSE of the model using equation (8). 

The contour and surface maps of the models and the 
known, as well as the gravimetric geoid heights of the 
points at various measurement point numbers, were 
plotted with surfer 11 to present their shape graphically. 

III. RESULTS PRESENTATION AND ANALYSIS 

Figures 3 and 4 respectively present the contour and 
surface maps of the third-degree model using 23 
observation points. Also, figures 5 and 6 respectively 
show the contour and surface maps of the known geoid 
heights of the 23 points. This was done to present 
graphically and compare the shapes of the model and the 
known geoid heights of the points to determine the 
resemblance as well as the variations between the model 
and the known geoid heights of the points. It can be 
respectively seen from figures 3 and 5 and figures 4 and 6 
that the contour and the surface maps of the third-degree 
model and the known geoid heights of the points are not 
identical which implies that the variations between the 
third degree model geoid heights and the known geoid 
heights of the points are considerably large. Thus, using 
23 points which is far larger than the third-degree surface 
10 terms, the accuracy of the model is very low. 



Fig. 4: Surface Map of Known 
Geoid Heights Using 23 Points 
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Fig. 8: Surface Map of Known 
Geoid Heights Using 22 Points 


Fig. 10: Surface Map of Third 
Degree Model at 22 Points 
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Fig. 6: Surface Map of Third 
Degree Model at 23 Point 


respectively seen from figures 7 and 9 and figures 8 and 
10 that the contour and the surface maps of the third- 
degree model and the known geoid heights of the points 
are not identical which also implies that the variations 
between the third degree model geoid heights and the 
known geoid heights of the points are very much large. 
Thus, using 22 points which is also far larger than the 
third-degree surface 10 terms, the accuracy of the model 
is also very low. 


Fig. 5: Contour Map of Third 
Degree Surface at 23 Points 

Figures 7 and 8 respectively present the contour and 
surface maps of the third-degree model using 22 
observation points. Also, figures 9 and 10 respectively 
show the contour and surface maps of the known geoid 
heights of the 22 points. This was also done to present 
graphically and compare the shapes of the model and the 
known geoid heights of the points to determine the 
resemblance as well as the variations between the model 
and the known geoid heights of the points. It can be 


Fig. 7: Contour Map of Known 
Geoid Heights Using 22 Points 


Fig. 9: Contour Map of Third 
Degree Surface at 22 Points 
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Again, figures 11 and 12 respectively show the contour 
and surface maps of the third-degree model using 21 
observation points. Also, figures 13 and 14 respectively 
present the contour and surface maps of the known geoid 
heights of the 21 points. This was as well done to present 
graphically and compare the shapes of the model and the 
known geoid heights of the points to determine the 
resemblance as well as the differences between the model 
and the known geoid heights of the points. It can be 



346000 348000 350000 352000 354000 356000 358000 360000 362000 

Fig. 11: Contour Map of Known 
Geoid Heights Using 21 Points 



Fig. 13: Contour Map of Third 
Degree Surface at 21 Points 


Figures 15 and 16 respectively present the contour and 
surface maps of the third-degree model using 11 
observation points. Also, figures 17 and 18 respectively 
show the contour and surface maps of the known geoid 
heights of the 11 points. This was as also done to present 
graphically and compare the shapes of the model and the 
known geoid heights of the points to determine the 
resemblance as well as the variations between the model 
and the known geoid heights of the points. It can as well 
be respectively seen from figures 15 and 17 and figures 
16 and 18 that the contour and the surface maps of the 


respectively seen from figures 11 and 13 and figures 12 
and 14 that the contour and the surface maps of the third- 
degree model and the known geoid heights of the points 
are not identical which implies that the differences 
between the third degree model geoid heights and the 
known geoid heights of the points are greatly large. Thus, 
using 21 points which are as well larger than the third- 
degree surface 10 terms, the accuracy of the model is well 
very low. 



Fig. 12: Surface Map of Known 
Geoid Heights Using 21 Points 



Fig. 14: Surface Map of Third 
Degree Model at 21 Points 


third-degree model and the known geoid heights of the 
points are identical which implies that the variations 
between the third degree model geoid heights and the 
known geoid heights of the points are very small. Thus, 
using 11 points which is very close to the third-degree 
surface 10 terms, the accuracy of the model is very high. 
This also implies that the accuracy of the local geometric 
geoid model is highest when the number of points used is 
either almost or equal to the number of polynomial 
geometric model terms. 
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Fig. 15: Contour Map of Known 
Geoid Heights Using 11 Points 



Fig. 17: Contour Map of Third 
Degree Surface at 11 Points 


Also, figures 19 and 20 respectively present the contour 
and surface maps of the fifth-degree model using 23 
observation points. Again, figures 21 and 22 respectively 
show the contour and surface maps of the known geoid 
heights of the 23 points. This was as well done to present 
graphically and compare the shapes of the model and the 
known geoid heights of the points to determine the 
resemblance as well as the variations between the model 
and the known geoid heights of the points. It can also be 
respectively seen from figures 19 and 21 and figures 20 
and 22 that the contour and the surface maps of the fifth- 



Fig. 16: Surface Map of Known 
Geoid Heights Using 11 Points 



Fig. 18: Surface Map of Third 
Degree Model at 11 Points 


degree model and the known geoid heights of the points 
are approximately identical which implies that the 
variations between the fifth-degree model geoid heights 
and the known geoid heights of the points are small. That 
is, using 23 points which is close to the fifth-degree 
surface 21 terms, the accuracy of the model is slightly 
higher. This also means that the accuracy of the local 
geometric geoid model increases as the number of 
measurement points approximates the number of 
geometric model terms. 



Fig. 19: Contour Map of Known 
Geoid Heights Using 23 Points 



Fig. 20: Surface Map of Known 
Geoid Heights Using 23 Points 
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Fig. 21: Contour Map of Fifth 
Degree Surface at 23 Points 


Besides, figures 23 and 24 respectively present the 
contour and surface maps of the fifth-degree model using 

22 observation points. Also, figures 25 and 26 
respectively show the contour and surface maps of the 
known geoid heights of the 22 points. This was also done 
to present graphically and compare the shapes of the 
model and the known geoid heights of the points to 
determine the resemblance as well as the differences 
between the model and the known geoid heights of the 
points. It can also be correspondingly seen from figures 

23 and 25 and figures 24 and 26 that the contour and the 



Fig. 23: Contour Map of Known 
Geoid Heights Using 22 Points 



Fig. 25: Contour Map of Fifth 
Degree Surface at 22 Points 



Fig. 22: Surface Map of Fifth 
Degree Model at 23 Points 


surface maps of the fifth-degree model and the known 
geoid heights of the points are very much identical which 
implies that the variations between the fifth degree model 
geoid heights and the known geoid heights of the points 
are truly small. That is to say, using 22 points which is 
very close to the fifth-degree surface 21 terms, the 
accuracy of the model is very high. This also implies that 
the accuracy of the local geometric geoid model tends to 
the highest as the number of measurement points is 
closest to the number of geometric model terms. 



Fig. 24: Surface Map of Known 
Geoid Heights Using 22 Points 



Fig. 26: Surface Map of Fifth 
Degree Model at 22 Points 
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As well, figures 27 and 28 respectively present the 
contour and surface maps of the fifth-degree model using 
21 observation points. Also, figures 29 and 30 
respectively show the contour and surface maps of the 
known geoid heights of the 21 points. This was over 
again done to present graphically and compare the shapes 
of the model and the known geoid heights of the points to 
determine the resemblance as well as the variations 
between the model and the known geoid heights of the 
points. It can again be respectively seen from figures 27 
and 29 and figures 28 and 30 that the contour and the 
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Fig. 27: Contour Map of Known 

Geoid Heights Using 21 Points 
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Fig. 29: Contour Map of Fifth 

Degree Surface at 21 Points 

Table 2 and figure 31 present the accuracy of the third- 
degree model at various measurement point numbers. 
This was done to compare the computed accuracy of the 
third-degree polynomial surface at various measurement 
point numbers. The smaller the computed Root Mean 
Square Error, RMSE, the better the accuracy of the 
model. It can be seen from table 2 and figure 31 that the 
accuracy of the third-degree geometric model at 23, 22, 
21 and 11 measurement points are correspondingly 


surface maps of the fifth-degree model and the known 
geoid heights of the points are extremely identical which 
implies that the variations between the fifth-degree model 
geoid heights and the known geoid heights of the points 
are extremely small. Therefore, using 21 points which is 
equal to the fifth-degree geometric model 21 terms, the 
accuracy of the model is highest. This also shows that the 
accuracy of the local geometric geoid model is highest 
when the number of measurement points is equal to the 
number of geometric geoid model terms. 



Fig. 28: Surface Map of Known 
Geoid Heights Using 21 Points 



Fig. 30: Surface Map of Fifth 
Degree Model at 21 Points 


0.632m, 0.6358m, 0.6358m and 0.034m. This means that 
the accuracy of the model is highest at 11 measurement 
points. This is as the number of points used (11) is very 
close to the number of the third-degree geometric model 
terms (10). This again, implies that the accuracy of the 
geometric geoid model tends to the highest as the number 
of measurement points approximates the number of the 
model terms. 
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Table 2: Accuracy of the Third Degree Model at Various 
Measurement Point Numbers 
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rd Degres-10 Terms Geoid Surface 


23 Points 

22 Points 

21 Points 

11 Points 

RMSE (m) 
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0.034 


■ RMSE (m) 
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Fig. 31 :Plot of Accuracy of the Third Degree Model 
at Various Measurement Point Numbers 

Again, table 3 and figure 32 present the accuracy of the 
fifth-degree model at various measurement point 
numbers. This was as well done to compare the obtained 
accuracy of the fifth-degree polynomial geoid model at 
various measurement point numbers. The smaller the 
computed RMSE, the better the accuracy of the model. It 
can be seen from table 3 and figure 32 that the accuracy 
of the fifth-degree geometric model at 23, 22 and 21 
measurement points are respectively 0.4333m, 0.0046m 
and 0.0003m. This means that the accuracy of the model 
is highest, less than 1mm at 21 measurement points. Also, 
at 22 observation points, the accuracy of the model is 
within 5mm. This once more implies that the accuracy of 
the polynomial geometric geoid model tends to the 
highest as the number of measurement points 
approximates the number of the model terms and in a 
unique solution, when the number of observation points is 
equal to the number of the polynomial geometric model 
terms, the model accuracy is highest. 


Table 3: Accuracy of the Fifth Degree Model at Various 
Measurement Point Numbers 
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Fig. 32: Plot of Accuracy of the Third Degree Model 
at Various Measurement Point Numbers 

IV. CONCLUSION AND RECOMMENDATIONS 

1. This paper has determined the relationship between 
polynomial geometric geoid surface terms and 
observation points number and effect in the accuracy 
of geometric geoid models. 

2. The paper has presented that the accuracy of the 
geometric geoid model tends to the highest as the 
number of measurement points approximates the 
number of the model terms. 

3. The paper has also presented that the accuracy of the 
polynomial geometric geoid model is highest when 
the number of observation points is equal to the 
number of the polynomial geometric model terms 

4. The obtained results have shown that the least 
squares model for fitting polynomial geometric 
surfaces to geoid heights should not just be the 
observation equal to the estimate but the residual 
should be considered as it can be used as a check. 
This is for the reason that the residual obtained from 
the difference between the estimate and the 
observation is equal to the difference between the 
model geoid height and its corresponding known 
geoid height. 

5. The paper has also recommended that the geometric 
method of local geoid model determination should 
be strictly applied in small areas. Where the method 
will be applied in considerable large areas, having 
considered the variation of the geoid heights, higher 
degrees polynomial geometric surfaces with larger 
numbers of terms approximating the number of 
observation points should be applied. This will 
enable the proper fit of the polynomial surface to the 
known geoid heights, as well as high accuracy to be 
obtained. 
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